Galaxy Evolution 


One can imagine at least 5 types of galaxy evolution: photo- 
metric evolution (i.e., changes in galactic luminosity, color, and 
spectra), chemical evolution, morphological evolution, dynamical 
evolution (i.e., changes in the types of stellar orbits), and num- 
ber evolution (i.e., changes in the number of galaxies present in 
the universe). The last three are difficult to study analytically, 
since they are closely related. For example, collisions between 
galaxies may turn spirals into ellipticals, and consume smaller 
galaxies. At the same time, new galaxies may form from debris 
thrown out in a collision. In addition, morphological evolution of 
spirals can be caused by environmental effects: in a rich cluster, 
ram-pressure sweeping of the interstellar medium may occur as 
the galaxy passes through x-ray gas. 





Fortunately, photometric and chemical evolution are more easily 
handled. Each can be parameterized in a fairly straightforward 
manner. 


Photometric Evolution 


[Tinsley 1980, Fund. Cosmic Phys., 5, 287] 

[Rana 1991, ARAA, 29, 129] 

[Renzini & Buzzoni 1986, in Spectral Evolution of Galaxies, p. 195] 
[Bruzual & Charlot 2003, MNRAS, 344, 1000] 

[Maraston 2005, MNRAS, 362, 799] 


The photometric evolution of galaxies is a particularly straight- 
forward problem. There are several generally available computer 
codes that can trace the photometric evolution of a simple (or 
complex) stellar population throughout cosmic history. However 
to first order, it is possible to estimate such properties analyti- 
cally. 





In order to make the mathematics a bit more straightforward, 
let’s choose a fiducial mass, mı, which can be the mass of the 
Sun. Such a star will have a main-sequence luminosity, 41, and 
a main-sequence lifetime of rı. We will then define the variable, 
m as the dimensionless mass of a star, i.e., for a star of mass M, 
m= M / mı. 








First we need is an Initial Mass Function (IMF) of stars, which 
gives the number of stars born as a function of mass. The IMF 
is usually written as a power law, t.e., 


o(m)dm = ġpım tt dm (8.01) 





where x describes the steepness of the function. [Most applica- 
tions use the local neighborhood IMF slope of x = 1.35 found by 
Salpeter (1956), but more and more people are going to broken 
power laws, such as those of Miller & Scalo (1979), Kroupa et al. 
(1993), or Chabrier 2003).] The constant ¢, is simply there for 
normalization purposes, so that a star cluster (or galaxy) with 
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total mass Mo satifies 


Mo = Mo [ mo(m) dm (8.02) 


The second piece of information we need is the mass-luminosity 
relation for main-sequence stars. This is well known from stellar 
evolution, and, although the slope of this relation does change 
along the main sequence, for our purposes, we can approximate 
the function as a single power law, with slope a. Thus, 


lg = lym (8.03) 


where the subscript d represents the luminosities of dwarfs (i.e., 
main sequence stars). While the mass-luminosity relation of the 
main sequence isn’t quite a single power law, the approximation 
isn’t too bad if you take a ~ 3.5. 


The third item needed is the length of time stars spend on the 
main sequence. The lifetime of a star is simply proportional to 
the energy available to the star divided by rate at which the star 
is emitting its energy, i.e., its luminosity. Since the available 
energy is proportional to the available mass, and the luminosity 
is proportional to mass (though the mass-luminosity relation), 
m m 
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Alternatively, this equation can be inverted. After a time t, the 
turnoff mass of a single-age stellar population will be 


one t 1/(1-a) 
Min X tl SO) Min = (=) (8.05) 
1 





The final pieces of information that are needed are the lifetime of 
a typical post-main sequence (giant) star (Tg), the average lumi- 
nosity of a giant star, 4, and the mass of a typical white dwarf, 
w. All of these can be derived (approximately) from models of 
stellar evolution. Note that Tg < Ta, so the stellar main-sequence 
lifetime is also the star’s total lifetime. 


We will now consider the evolution of a set of stars all born (with 
the same metallicity) at the same time. This is called a Simple 
Stellar Population (SSP). The properties of more complicated 
systems can be inferred from the summation of several SSPs. 


Luminosity Evolution 


Let’s first calculate the luminosity evolution of a cluster (or 
galaxy) of stars with total mass Mo all born at the same time. 
The total luminosity of main sequence stars is easy to compute: 
all we have to do is sum up the luminosity of all stars on the 
main sequence. The lower end of the sum is the minimum mass 
of an energy-generating star; the upper end is defined by the 
main sequence turnoff. In other words, 
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Note that since the exponent is negative, and low mass stars 
live (essentially) forever, the last term in the above equation is 
negligible. Thus, the total luminosity of dwarf stars is 


Lp(t) 
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Calculating the total luminosity of giant stars is equally simple. 
First, note that the timescale for giant branch evolution is much 
faster than that for main sequence evolution. Thus, the key is 
to estimate the number of stars currently turning into giants; 
when this number is multiplied by the length of time a typical 
star remains a giant, and the mean luminosity of the star, the 
result is the total giant star luminosity. Now consider: the rate 
at which main sequence stars turn into giants is defined by how 
many stars are at the main sequence turnoff, and how much of the 
main sequence is eaten away per unit time interval. According 
to our definition of the initial mass function, the number of stars 
at the main-sequence turnoff is 
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and the number of stars turning into red giants during a time At 
is 
_ aN dN dm 
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so the total luminosity of giant stars is 
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If we now substitute for t for m using (8.05) take the derivative, 
and combine terms, the result is 
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Thus, the total luminosity of the stellar system, as a function of 
time, is 
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This can be simplified a bit if we define G(t) as the ratio of giant 
star luminosity to dwarf star luminosity, t.e., 


GG 226 a E (Ly (8.13) 


ee (a-lan n 


Since a > 1, the exponent of time in (8.13) is significantly less 
than one. Hence writing the expression for G(t) in this way 
illustrates that the variable is only a weak function of time. This 
expression can then be further simplied by using (8.03) and (8.05) 
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Note that the first part of this equation is very close to one. 
Furthermore, the term in brackets is simply the ratio of the total 
energy generated by the star on the giant branch to the total 
energy the star generated on the main sequence. Since stars turn 
off the main sequence when ~ 10% of their total fuel is gone, 
but eventually burn ~ 70% of their fuel, G(t) ~ 6. (From an 
observers point of view, it’s a bit more complicated, since most 
of the light a giant star produces is red, while the main sequence 





light may be blue. Thus, the exact value of G(t) depends on 
the wavelength of observation: in the blue, G(t) ~ 1, but in the 
red, G(t) > 10.) Using this notation, the luminosity of a stellar 
population, as a function of time, is 
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Cluster Evolution and qo 


A traditional method of measuring go, which dates back (in idea) 
all the way to Hubble, is to use the brightest elliptical galaxies as 
a standard candle, and look for curvature in the Hubble diagram. 
Recall from the very first homework assignment that, when the 
equation for distance modulus in an Einstein-de Sitter universe 
is expanded out in a power series, 


Qo 
Mpol = Mboi + 5 log + 1.086(1 = =) z (8.17) 
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So, as z increases, the relation between Mpo; and log z increas- 
ingly departs from a straight line Hubble law, in a manner that 
depends on qo (or Qo and A). In the A = 0 case, high-z galaxies 
will appear brighter than they should (according to the linear 
Hubble law) by 


Q 
Am = 1.086(1 — a (8.18) 


Now consider what passive stellar evolution will do to the Hubble 
diagram. Let’s assume that all the stars in an elliptical galaxy 
were born (essentially) at the same time. From (8.13) 
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We can find the derivative of G with respect to t via (8.13). If we 
do that, then do a little mathematical manipulation, we arrive 


at 
dln £ 1 G 
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With a Salpeter mass function, a ~ 3.5, and G = 6, E x —0.5. 


How bad is this? For an Einstein-de Sitter universe, the lookback 
time for this galaxy is 
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where we have simply expanded the equation for lookback time 
in a power series near z = 0. Using this, the luminosity evolution 
that has occurred between redshift z and today is 
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If we convert ln £ to magnitude, and again expand the resulting 
equation in a Taylor series about z = 0 (i.e., using ln(1 +x) ~ x), 
then 
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The above equation says that at large z, galaxies will appear 
much brighter than they should according to the Hubble law. 
Moreover, if you set (8.18) and (8.24) equal 


Q Q 
Am = 1.086(1- -> z) = 1.086Ez => a =1—-E ~ 1.5 (8.25) 


In other words, if you want to use elliptical galaxies as standard 


candles, you must correct for the passive evolution of their stars, 
and the amount of this correction is large! 


AQ 9 = Qo(observed) — Qo(true) = 3 (8.26) 
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Photometric Mass-to-Light Ratio 


In addition to a population’s total luminosity, there are several 
other quantities which can be computed from first principles. 
First is the photometric mass-to-light ratio. Recall the total lu- 
minosity of main sequence stars in a stellar population is 





oe Moos f {mp * —mz-*} (8.06) 
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where min is the turnoff mass of the main sequence. As stated 

previously, a — x is usually greater than zero, so the last term of 

this equation is negligible; stars at the low end of the mass se- 

quence do not contribute significantly to the system’s luminosity. 


On the other hand, low mass stars can be an important part of 
the cluster’s mass. From (8.02), the visible mass of stars in a 
galaxy is given by 
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(neglecting the mass in invisible stellar remnants). If « > 1 
(as is normally assumed), the exponents in (8.27) are less than 
zero, and the value of M is dominated by the last term in the 
equation. Physically, this means that most of the stellar mass of 
a population resides in low mass stars. 


The combination of (8.06) and (8.27) means that it is virtually 
impossible to determine a population’s photometric mass-to-light 
ratio from observations alone. One can always drive up M/L by 
postulating the existence of low-mass stars which add mass (via 
8.27), without adding luminosity (8.06). 
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Stellar Evolutionary Flux 


An interesting property associated with stellar populations is the 
stellar evolutionary flux, i.e., the number of stars passing through 
any (post main-sequence) phase of evolution at a given time. Just 
as the amount of water flowing down a river is controlled by the 
spillway of a dam, the number of stars evolving through various 
stages of evolution is controlled by the rate at which these stars 
move off the main sequence. (All other rates are much faster 
than this.) As we have seen previously, this rate is 
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From (8.05) 
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so the stellar evolutionary flux is 
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Naturally, the number of stars evolving through any phase of 
evolution in a galaxy is proposal to the number of stars in the 
galaxy. The best way to take this dependence out is to normal- 
ize the evolutionary flux to the size of the galaxy. This cannot 
be done with mass, since we do not know the total number of 
stars that are present in the system. However, we do know the 
population’s luminosity, which is defined in (8.16). So, when we 
normalize (8.30) to this luminosity, we arrive at the population’s 
luminosity-specific stellar evolutionary flux. After a bit of 
math, this quantity comes out to be 
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Note that the exponent of time is less than 1; thus B is rather 
insensitive to the age of the stellar population. (The difference 
between the stellar evolutionary flux of a 7 Gyr stellar population 
and a 12 Gyr population is less than a 25%.) Similarly, the stellar 
evolutionary flux is not very sensitive to the exact value of z. 








Detailed calculations show that the luminosity-specific stellar 
evolutionary flux is B ~ 2 x 107"! stars yr7! De. for all old (or 
moderately old) stellar populations. This makes is easy to pre- 
dict the number of any post-main sequence type of star present 
in a galaxy: since you know B, and you can measure the lu- 
minosity of the galaxy, all you need is the lifetime of the phase 
in question. For example, the lifetime of a planetary nebula is 
T ~ 25,000 years. If a galaxy is observed to have an absolute 
luminosity of L = 10'4Lo, then the number of PN in the galaxy 
is N(PN) ~ B - £ - T = 50, 000. 





Note: a common idea these days is that all (or many) planetary 
nebulae are formed out of binary star interactions. If so, then 
the previous calculation may not work, since it does not take into 
account binary-star interactions. (Interestingly, however, the es- 
timate does appear to come close to predicting PN numbers.) 
Similarly, if not all stars evolve through a given stage, then the 
calculation may be in error. 
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The plot above shows how a population’s bolometric-luminosity 
specific stellar evolutionary flux changes with the systems age 
and IMF. Note that the x-axis (population age) is a log quantity, 
while the y-axis is linear. All old (> 3 Gyr) stellar population 
have a rate of ~ 2 x 10!! stars per year per solar luminosity, 
independent of age, IMF, or metallicity. 


[Renzini & Buzzoni 1986, in Spectral Evolution of Galaxies, p. 195] 
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Mass Loss from Stars 


Another useful quantity to know is the rate of mass loss from 
stars as a function of time. As before, the key to calculating 
this is to realize that almost all the mass lost from stars comes 
during the post main-sequence phase of evolution. Thus, the rate 
at which a population of stars loses mass is just by the rate at 
which stars move off the main sequence times the amount of mass 
each star loses. If m is the initial mass of the star, and w the 
mass of the stellar remnant (i.e., white dwarf), then the mass 
ejection rate is 








E(t) = Nin(min — w) = My (Min) Z (men — w) (8.32) 


After performing the same substitutions as was done for the cal- 
culation of stellar evolutionary flux, this simplifies to 
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Once again, if we normalize to luminosity via (8.16), we can 
derive the luminosity specific mass loss rate 
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Numerically, this works out to a mass loss rate of ~ 0.02Mo 
per Gigayear per unit solar luminosity for an old (10! yr) stellar 
population. When integrated over the lifetime of a galaxy, the 
result is that ~ 15% of the initial stellar mass will be lost during 
a Hubble time. 
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More Complex Systems 


Real galaxies usually contain multiple stellar populations, each 
with its own mass, age, and metallicity. Computing the behavior 
of such a system simply means summing variations SSP compo- 
nents. It would, however, be useful to some zeroth order approx- 
imation to know the mix of stellar populations within a galaxy. 





A common way to do this is to simply assume that the star- 
formation rate of a galaxy declines exponentially with time, with 
some decay rate T, t.e., 


M(t) « e7*/7 (8.35) 


Irregular and late-type spiral galaxies are assigned large values 
of T, suggesting an almost constant star-formation rate over the 
history of the universe. Elliptical galaxies have small values of 
T, thus reflecting the fact that their star-formation has e-folded 
away long ago. Note however, that, many/most galaxies likely 
undergo a succession of discrete starbursts, with low-level star 
formation occuring in between times of extreme activity. So these 
models may be of only moderate use. 
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Problems in Photometric Evolution Calculations 


The above analytic solutions are only useful for guidance. To 
compare with observations, predictions must be made in specific 
bandpasses, or for specific absorption lines. This requires numer- 
ical calculations which include 


sets of stellar isochrones, detailing the precise number of stars at 
any position in the HR diagram as a function of age. 


sets of model stellar atmospheres, giving the emitted spectrum 
(or broadband color) of each star, as a function of its tempera- 
ture, surface gravity, and metallicity. 


an understanding of various non-stellar processes that may effect 
the emergent spectrum from galaxies. Such effects include inter- 
nal extinction due to dust, the presence of emission-lines from 
ionized gas, and (for high-redshift objects) absorption due to the 
intervening intergalactic medium (i.e., the Lya forest). 





Today, the most controversial parts of these codes include their 
treatment of thermally pulsing AGB stars, and the proper mod- 
elling of the horizontal branch and post-AGB stars. The former is 
important since it involves extremely luminous stars, which, due 
to their mass loss, may be obscured by circumstellar extinction. 
The latter is important for modeling old stellar systems which 
have very few blue stars. For these objects, a small change in 
the morphology of the horizontal branch (i.e., the ratio of red 
to blue objects), or the number of slow-evolving post-AGB stars 
can change the system’s UV color dramatically. 
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